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Abstract. We demonstrate a close relation between Coulomb effects in non-local electron 
^ I transport and non-local shot noise in three-terminal metallic conductors. Provided the whole 

structure is normal, cross-correlations in shot noise are negative and Coulomb interaction tends 
to suppress both local and non-local conductances of the system. The behavior of normal- 
superconducting- normal structures at subgap energies is entirely different. In the tunneling 
2 ' limit non-local differential conductance of such systems are found to have an S-like shape and 

can turn negative at non-zero bias. At high transmissions crossed Andreev reflection yields 
positive noise cross-correlations and Coulomb anti-blockade of non-local electron transport. 

o 

1. Introduction 

^1 It is well known that discreteness of electron charge stays behind a number of fundamentally 

f^ ■ important physical phenomena such as, e.g., shot noise in mesoscopic conductors [Ij and 

OO I Coulomb blockade of charge transfer in tunnel junctions [2]. Exactly for this reason these two 

seemingly different phenomena turn out to be closely related to each other: Coulomb blockade is 
stronger in conductors with bigger shot noise [31 H] . This fundamental relation was subsequently 
^— s I confirmed in experiments [5j. Later a close relation between shot noise and Coulomb blockade 

^SJ ■ was also identified in hybrid normal-superconducting (NS) structures [6], where doubling of 

elementary charge due to Andreev reflection becomes important at subgap energies. 

Recently the same idea was extended [7] to non-local effects in electron transport across three- 
terminal normal-superconducting-normal (NSN) systems where entanglement between electrons 
in different normal terminals can be realized. Non-local electron transport in such systems is 
^.' determined by an interplay between elastic cotunneling (EC) and crossed Andreev reflection 

(CAR) and was recently investigated both experimentally ^ [9l [lOl [TTl [12] and theoretically 
|131 \T^ 115] , While non-interacting theory predicts that CAR never dominates over direct 
electron transfer, both positive and negative non-local signals have been detected in a number 
of experiments [SI [9l [Til [12] • Theoretically it was argued that CAR could prevail over EC 
in the presence of Coulomb interactions [16j or an external ac field [17]. Negative non-local 
conductance was also predicted in interacting single-level quantum dots in-between normal and 
superconducting terminals [15] . 

In Ref. [7j we have already demonstrated that interaction effects in non-local transport and 
non-local shot noise in NSN systems are intimately related. This relation, however, turns out 
to be much more complicated than in the local case [31 [H [S|- The main reason for that is 
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Figure 1. Schematics of the system mider consideration. 

positive cross-correlations in shot noise which may occm' in normal-superconducting hybrids 
[U [19] . In NSN structures such positive noise cross-correlations were demonstrated theoretically 
[201 [211 [7] and experimentally [22] . Note that this feature is specific to superconducting systems 
and is totally absent in normal ones where cross-correlations of fluctuating currents are known 
to always be negative |lj. Hence, it would be interesting to extend our theory of non-local 
electron transport in NSN systems in the presence of electron-electron interactions [7] to normal 
conductors and compare the corresponding results derived for superconducting and normal 
structures. This is the main goal of the present paper. 

Our paper is organized as follows. In Sec. 2 we define our model and outline the key steps 
of the derivation of the effective action for our three-terminal metallic structure. In Sec. 3 we 
re- formulate our results in terms of equivalent Langevin equations describing real time dynamics 
of fluctuating voltages and currents and demonstrate how the expressions for both local and non- 
local conductances in the presence of Coulomb interaction can be related to the corresponding 
shot noise correlators. Sec. 4 and 5 are devoted to the analysis of the effect of electron-electron 
interactions on the conductance matrix respectively for superconducting and normal central 
electrodes. A brief summary of our main observations is presented in Sec. 6. 

2. The model and effective action 

Let us consider a thin wire, which can be either normal or superconducting and is connected 
to two normal metallic leads via two barriers. Quite generally, these barriers are characterized 
by two sets of channel transmissions T^ , T^ . For the sake of simplicity here we assume that 
the first junction has A^i conducting channels with the same transmission Ti, while the second 
one — A'^2 channels with transmission T2. Provided the wire is in the superconducting state, its 
energy spectrum has the gap A. This gap, of course, equals to zero should the wire be in the 
normal state. In the latter case it is convenient to define the normal state conductances and the 
Fano factors of the barriers as follows 
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If the wire is superconducting, we define Andreev conductances and Andreev Fano factors, which 
determine the transport properties across NS barriers in the subgap regime eVr,T <^ A: 
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where r^ = T^/{2 — Tj.)'^ represent effective Andreev transmissions of NS barriers. Both barriers 
are supposed to have capacitances Ci, C2 and the leads 1 and 2 are characterized by large Ohmic 



conductances Gf^ and Gl*^. We also define the corresponding dimensionless conductances of the 
electromagnetic environment 
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The corresponding dimensionless conductances in case of a normal wire g^ may differ from g^ 
due to an additional contribution from the wire resistance. 

Weak electromagnetic coupling between two NS barriers (e.g. via modes propagating in 
the superconductor |16) ) will be disregarded. Our main goal here is to evaluate electric currents 
A(^i) ^) and /2(^i5 ^) across the barriers 1 and 2. Below we will consider two most interesting 
limits: (i) the subgap regime T, eVi, eV^ <C A in the case of a superconducting wire, and (ii) the 
regime of high energies or a normal wire. 

The general Hamiltonian of our system can be expressed in the form [7] 

H = Hi+H2 + //wire + Ht,1 + Ht,2, (4) 

where 

Hr= J2 J ^^V^r,a ( "^ " /^ ) '^r,a, T = 1,2, 

are the Hamiltonians of the normal metals, m is electron mass, ^ is the chemical potential, 

//wire = dx J2^^o,{-^- f^jXa + ^x\xl + ^*)aXt 

is the Hamiltonian of a superconducting wire and 

are tunneling Hamiltonians describing transfer of electrons across the contacts with area Ar and 
tunneling amplitude tr- As we already indicated above, we will assume that both barriers 
are uniform implying that all A^^ = kpAr/'^TT conducting channels in the r-th barrier are 
characterized by equal transmission values 

2 



Tr = 47r^MrZ^wire|irP/(l + VT^ i^r l^wire | ^r P j , (6) 



where Vj (j = 1,2, wire) is the density of states in the corresponding terminal. Finally, we 
note that fluctuating phases fr introduced in Eq. Q are linked to the voltage drops across the 
barriers Vr by means of the standard relation (pr = eVr and are treated as quantum operators. 

We will proceed in a standard manner and eliminate fermionic variables expressing the kernel 
J of the Keldysh evolution operator via path integral over the phase fields [21 [23] 

J= f Yl Vip^V^^ expCzSenviv'f , ^'^] + iStI^^., ^r]). (7) 

■^ r=l,2 

where ^p^ and ^p^ are fluctuating phases defined respectively on the forward and backward 
branches of the Keldysh contour, S'env is the action of electromagnetic environment and the 
term iSt accounts for electron transfer between the terminals. In the case of linear Ohmic 
environment considered here one has [2] 

T-=l,2 



where ipr = {<f^ + ipf)/2, ip^. = ip^ — ip^ and 
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The term iSx is derived from the tunnel Hamiltonians ([5]) and reads 

/ G^' h \ 
iST = iT\nG-\ Q^=\ A G^re h , (9) 

V 4 G^' J 

where 4x4 matrices G~ represent the inverse Keldysh Green functions of isolated normal leads 
(j = 1, 2) and the wire, and tr is diagonal 4x4 matrix in the Nambu - Keldysh space 



t 



r 



I -tre-''^r 














tre-'^fi? 














Ue'^Pr 





l 








-Ue'^? 



(10) 



After some exact manipulations we obtain 

iSt = tr In [i - tlGiiiGwire - ilGihG^ire] ■ (11) 

While the expression (jlip for the action remains formally exact it is still too complicated to 
be directly employed in our calculations. In order to proceed we will make several additional 
steps which yield significant simplifications. First, we restrict ourselves to the limit of high 
conductances 

5^, g^ » 1, (12) 

in which case phase fluctuations are weak and it suffices to expand the action pip to the second 
order in ip~ , cf., e.g., [31E1IM]- Next we assume that the resistance of the wire segment between 
the junctions, ri2, is small as compared to the junction resistances, ri2 ^ 1/G^. In this case 
one can additionally expand the action in powers of the wire Green function connecting the 
two junctions. In the case of a superconducting wire and at T, eVr <^ A the corresponding 
calculation was elaborated in Ref. [7]. In this case the effective action of our system can be cast 
to the form 

iSr = iSii + ^5*22 + iSu, (13) 

where the contributions iSn and iSi2 read 

^Sn = -^^ J dt^.p^^ - J dtdt' ''^^'^^^/^^''\ (14) 

,5i2 = i^Jdt {^W, + Vi2^r) - / dtdt'^^^^^^l-^^^ (15) 

while the term 1822 is obtained by interchanging the indices 1 o 2 in Eq. (|14p . The functions 
5** are defined as 

5**' = Gfi^M(t - t'){l - /3i + /3i cos[2(^f ]) + 2G^iM{t - t')(ai - m cos[2(y9f ]) 

+(Gf2^/2)M(t - t'){4 cos[<^f + <^«'] + /.r cos[(^f - (^f']), (16) 

5**2' = -G^iM{t - t'){l - /3i + /3i cos[2^f ]) - Gf2^M(t - t'){l -h + h cospc/;*/]) 

+ (Gf2^/2)M(t - t')(7+ cos[v9f + ^f] - 7- cos[<^f - <^«']). (17) 



Here we denoted 93** = ipr{t) — (pr{t')- Other parameters entering in Eqs. (fT6]) and p!7|) read 

ar = Tr{l - 2Tr)/^TiT2, r]r = 2Tr {I - Tr) / ^/ttT^ , 

4 = ±(4r, - 3) + 1/Vn^ {r = 1, 2), ^^ = ±1 + (1 - 2ri - 2t2 + 4riT2)/Vn^, (18) 

while zero bias non-local conductance has the form 

r<NSr<NS 
Gf/ = ^iV^^I)(2.A,n,r2). (19) 

Here D{u}, r, r') is the diffuson, which is defined as a solution of the following diffusion equation 

-iLo-DVl)D{u;,r,r') = d{r-r'). (20) 



In a simple quasi-one-dimensional geometry of Fig. 1 one finds [15] Gj^ = Gii G22 R^ e^ ' ^/2, 
where d is the distance between two NS barriers and R^ is the resistance of the piece of the 
wire with the length equal to the superconducting coherence length ^ = \JD j ^. It is important 
to emphasize that all order terms in tr are fully accounted for in Eqs. (|14p - (jl8p . i.e. the action 
applies for arbitrary transmission values T\^2 (or t\^2) ranging from zero to one. 

The calculation in the normal case goes along the same lines, however the resulting effective 
action turns out to be somewhat more complicated because of retardation effects related to 
diffusion of electrons between the barriers. For the sake of simplicity here we avoid presenting 
an explicit form of the effective action in the normal case. Rather we will proceed directly to 
the final expressions for non-local currents which follow from this action. The corresponding 
expressions will be presented in Sec. 5. 

3. Langevin equations and interaction correction to the current 

It is well known that the quadratic in (/3J~2 effective action can be exactly rewritten in terms 
of the corresponding Langevin equations |241 [25l [26] which describe the current balance in our 
system. In case of a superconducting wire these equations read [7] 

Cii;i + (Gf + Gfi^)z;i - Gf/z;2 = GfV^ + ^i + 6, 

C2V2 + iGf + G^i)v2 - Gf/t;i = GfV2 + ^f + ^2. (21) 

Here Vr are the fluctuating volatge drops across the junctions, ^f are stochastic variables with 
pair correlators 

(ef (i)C^*')) = GfM(t - t') (22) 

describing Gaussian current noise in the shunt resistors, while the variables ^^ with the 
correlators 

iUtMt')) = S',\' (23) 

describe shot noise in NS barriers. 

In case of normal wire the Langevin equation looks similar. One needs to replace G^ by G^. 
Besides that, the non-local conductance exhibits retardation effects, i.e. one should replace, for 
example, G12V2 — )• ^-^dt'Gi2{t — t')v2{t'). Finally, the correlator of the noises also differs from 
(|23p . As we already pointed out, the corresponding expressions are rather cumbersome and for 
this reason are not presented here. 

Let us evaluate the current Ji across the first barrier. Solving Eqs. (|2ip perturbatively in 
l/S'r ^ 1, in the lowest non-trivial order in this parameter we get 

/i = Gfi^yi-Gf/y2-(6). (24) 



Here the average (^i) does not vanish since according to Eqs. (fT6|) . (fT7|) the noise ^i depends 
on the phases 991,2 , which, in turn depend on ^1,2 by virtue of Eqs. (|2ip . Hence, we obtain 



where the phase fluctuations 6ipr are determined from Eqs. (j2ip . We obtain 



(25) 



6Mt) = e I dt' —^ Cr (t') . 



(26) 



Here we assumed G12 ^ G-rr ^ G^'' and introduced the i?C— time tjic = Cr/Gf^. Substituting 
this expression into Eq. (j25p we get 
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Below we will make use of the above expressions and directly evaluate the non-local currents 
in both interesting limits of superconducting and normal wires. 

4. Superconducting wire 

We begin with a superconducting wire. In this case employing Eqs. ()16p . ()17p and performing 
the time integral in Eq. ()27p we arrive at the following expression for the current through the 
first barrier 



GnV, - G,2V2 - ^^li^l^i^l^Fo(2yi) + ^%|^Fo(2F2) 
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Taking the derivatives of the current (j28p with respect to Vi and V2 we define respectively local 
and non-local differential conductances of the first barrier. They read 
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(32) 



and ^(x) is the digamma function. The conductances (pO]) . ([3T]) are displayed in Fig. [2l 

We observe that both these conductances are affected by Coulomb interaction which yields 
non-trivial corrections to the corresponding non-interacting expressions. In the interaction 
correction to the local conductance in Eq. ()30p we recover the Coulomb blockade term [6] tx /3i 
and, in addition, three non-local contributions. The first of them oc r/i enhances the conductance, 
while the second one ex k^,7+ provides additional Coulomb suppression of dli/dVi. The last 
term oc K]^,7_ can be both positive (at ti^2 "^ 1) and negative (at bigger ri^2) implying the 
tendency to Coulomb anti-blockade in the latter case. The first term tx /32 in Eq. ()3ip for 
the non-local conductance has an opposite sign as compared to G12 (thus implying Coulomb 
blockade), while the second one oc Kf,^+ yields Coulomb anti-blockade. Finally, the third 
term oc kj~,7_ tends to suppress or enhance the absolute value of the non-local conductance 
respectively for k^ /qi^ + j-/g2^ > and k^ /g^^ + ■~f-/g2^ < 0. 

As we already argued [7] these non-trivial features of different Coulomb corrections are 
directly related to the corresponding contributions to shot noise. In particular, negatively 
and positively correlated noise terms are associated respectively with Coulomb blockade and 
anti-blockade terms in the above expressions for both local and non-local conductances. 
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Figure 2. Local (a) and non-local (b) differential conductances of a sample with 
superconducting wire, which are defined, respectively, in Eqs. ()30p and ()3ip . The parameters 
of the system are: T = 20 mK, G^^^ = 1 mS, G^^ = 3.2 mS, G^i^ = 10 fiS, g^^ = g^^ = 516, 
trc = W~^^ s. The transmissions of the junctions are small, Ti = 0.063, T2 = 0.11, the 
corresponding numbers of channels are, on the contrary, large, A'^i = A'^2 = 6.1 x 10^. 



5. Normal wire 

Let us now turn to the case of a normal metallic wire. In this case the current /i takes the form 
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where the normal state non-local conductance is defined as follows 
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Accordingly the differential resistances read 
dh ^N 2Gfi/3f 
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Here the functions Fi^iy) and -F(F) are again defined in Eqs. ([29]) and ([32 
The conductances (f35|) and ([36]) are depicted in Fig. [3l 
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Figure 3. Local (a) and non-local (b) differential conductances of a sample with normal wire, 
which are defined, respectively, in Eqs. ()35p and ()36p . The system parameters are the same as 
in Fig. El namely T = 20 mK, q^ = q^ = 516, T^ = 0.063, T2 = 0.11, iVi = iV2 = 6.1 x 10^, 
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mS, Gfs^ = 4.83 /uS. 



10 ^^ s. Accordingly, the normal conductances take the values Gii = 29.8 mS, G22 = 52 



Comparing these results with those obtained above for NSN systems we observe striking 
differences between them. The main difference is due to the fact that correlations in non-local 
shot noise in normal multi-terminal systems are always negative, while positive cross-correlations 
may occur in structures involving superconductors. Accordingly, in the normal case Coulomb 
interaction tends to always suppress non-local conductance, and the corresponding interaction 
correction in Eq. (p6|) depends on the voltage difference Vi — V2. In contrast, the expression for 
the non-local conductance in three-terminal NSN systems (|3ip contains Coulomb terms which 
depend both on Vi — V2 and Vi + V2 originating respectively from negative and positive cross- 
correlations in shot noise. As we already discussed, these terms describe respectively Coulomb 
blockade and Coulomb anti-blockade of the non-local conductance in the superconducting case. 



6. Discussion and summary 

In this paper we developed a theory elucidating a non-trivial physical relation between shot 
noise and Coulomb effects in non-local electron transport in three-terminal metallic structures. 
We evaluated non-local current-current correlators in such systems at arbitrary interface 
transmissions and arbitrary frequencies and directly related them to Coulomb effects in non-local 
electron transport. 

One of the important features of NSN systems under consideration is that in the tunneling 
limit almost no effect of Coulomb interactions on non-local conductance is expected if one of 



the applied voltages, Vi or V2, equals to zero. This effect is directly related to the cancellation 
between EC and CAR contributions to shot noise in the corresponding limit [20]. For nonzero Vi 
and V2 no such cancellation exists anymore and the non-local conductance dIi/dV2 approaches 
the S-like shape being enhanced at Vi ^ V2 and partially suppressed at Vi ^ —V2, see Fig. [2}3. 
Both these features have a clear physical interpretation. Indeed, at Vi ss — V2 negative cross- 
correlations due to EC dominate non-local shot noise leading to Coulomb blockade of non-local 
conductance while at Vi ~ V2 positive cross-correlations due to CAR prevail and Coulomb anti- 
blockade of non-local transport is observed. At higher interface transmissions only Coulomb 
anti-blockade of non-local conductance remains, which is again related to CAR-induced positive 
cross-correlations in shot noise. 

Comparing these results with those obtained in systems with normal central electrodes we 
observe striking differences. In particular, both local and non-local conductances always tend to 
be suppressed by Coulomb interaction and anti-blockade effects never occur in normal structures. 
The non-local Coulomb corrections to conductances depend on the voltage difference Vi — V2 
not on their sum ^1 + ^2, unlike in the superconducting case. All these features are directly 
related to the observation that only negative cross-correlations in shot noise occur in normal 
multi-terminal conductors [T]. 

It is interesting to point out that S-like shaped non-local signal predicted here was indeed 
observed in experiments with NSN systems |121 127] . A good agreement between our theory 
and the results [27] argues in favor of electron-electron interactions as a physical reason for 
the observed feature. Some of the features similar to those predicted here have also been 
observed in experiments [22]. It would be interesting to perform more experiments in systems 
under consideration both in superconducting and normal states and compare the results with 
our theoretical predictions. Extending experimental investigations to the normal case would 
hopefully allow for better characterization of the system parameters as well as for clearer 
demonstration of qualitative differences between normal and superconducting structures outlined 
above. 

Finally, we would like to make one more remark. In some cases non-linearities in both 
local and non-local differential conductances caused by Coulomb interactions may be combined 
with the zero bias anomalies resulting from the proximity-enhanced electron interference in 
diffusive normal leads [281 \29\ \30\ [15] . In this paper we disregarded this effect for the sake of 
simplicity. In practice it implies that here we considered the system with weakly disordered or 
sufficiently thick normal leads and sufficiently resistive barriers. If needed, zero-bias anomaly 
effects |28[ \29\ [301 1^5] can be included into our analysis in a straightforward manner. 
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